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ABSTRACT 

The issue of non-Gaussianity is not only related to distinguishing the theories of the 
origin of primordial fluctuations, but also crucial for the determination of cosmolog- 
ical parameters in the framework of inflation paradigm. We present a method for 
testing non-Gaussianity on the whole-sky CMB anisotropics. This method is based 
on the Kuiper's statistic to probe the two-dimensional uniformity on a periodic map- 
ping square associating phases: return mapping of phases of the derived CMB (similar 
to auto correlation) and cross correlations between phases of the derived CMB and 
foregrounds. Since phases reflect morphology, detection of cross correlation of phases 
signifies the contamination of foreground signals in the derived CMB map. The ad- 
vantage of this method is that one can cross check the auto and cross correlation of 
phases of the derived maps and foregrounds, and mark off those multipolcs in which 
the non-Gaussianity results from the foreground contaminations. We apply this statis- 
tic on the derived signals from the 1-year WMAP data. The auto-correlations of phases 
from the ILC map shows the significance above 95% CL against the random phase 
hypothesis on 17 spherical harmonic multipoles, among which some have pronounced 
cross correlations with the foreground maps. We find that most of the non-Gaussianity 
found in the derived maps are from foreground contaminations. With this method we 
are better equipped to approach the issue of non-Gaussianity of primordial origin for 
the upcoming Planck mission. 

Key words: cosmology: cosmic microwave background - observations - methods: 
data analysis 



1 INTRODUCTION 

The issue of non-Gaussianity in the cosmic microwave back- 
ground (CMB) has touched the most fundamental base in 
cosmology. It was first brought to attention by Ferreira, 
Magueijo and Gorski (1998) that non-Gaussian signal is 
present in the COBE data. Although it is almost certain 
that the departure from Gaussianity is induced by sys- 
tematic error (Banday, Zaroubi & Gorski 2000; Magueijo 
& Medeiros 2004), the discussions and focus about non- 
Gaussianity since then have been focusing primarily on pri- 
mordial origin. Mechanisms other than the simplest inflation 
model that have been proposed for primordial density fluc- 
tuations produce non-Gaussian fields (see Bartolo et al. 2004 
and references therein). Due to this reason, the issue of non- 
Gaussianity seems to be discussed separately from that of 
the determination of cosmological parameters. It is in the 
framework of inflation paradigm that the cosmological pa- 
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rameters can only be determined correctly from the angular 
power spectrum if the CMB temperature anisotropies con- 
stitute a Gaussian random field (GRF). Therefore, the issue 
of non-Gaussianity is not beyond the power spectrum, but 
still within the power spectrum. 

The statistical characterization of temperature fluctua- 
tions of CMB radiation on a sphere can be expressed as a 
sum over spherical harmonics: 

oo t 

AT(6,(p) = J2 ai m Yem(6,tp), (1) 

£ = 1 m=~i 

where ag m — \a^ m \ exp(i0£ m ). The strict definition of a ho- 
mogeneous and isotropic GRF, as a result of the inflation 
paradigm, requires that the moduli \ae m \ are Rayleigh dis- 
tributed and the phases <t>i m are uniformly random on the 
interval [0,27r]. The central limit theorem, however, guar- 
antees that a superposition of a large number of harmonic 
modes will be close to a Gaussian as long as the phases 
are random. Hence the random-phase hypothesis on its own 
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serves as a definition of Gaussianity (Bardeen et al. 1986; 
Bond & Efstathiou 1987). 

One of the most useful properties of GRF is that the 
second-order statistics, the 2-point correlation function or 
the angular power spectrum Ce 

(%m<y m '} = Cl <W<Lm' (2) 

furnish a complete description of the GRF. It is based on 
this analytically-simple but important property that the cos- 
mological parameters can be correctly determined from Ce- 
Accordingly, if non- Gaussian signals are present, either with 
primordial origin, or induced from data processing or sys- 
tematic error, the cosmological parameters derived from the 
Ce of such a "contaminated" field will have larger error bars. 
The issue of non-Gaussianity in CMB is therefore not only 
related to discriminating the theories of origin of primordial 
fluctuations, it is also fundamental, in the framework of in- 
flation paradigm, for the determination of the cosmological 
parameters. 

To test non-Gaussianity, the next order statistics: the 3- 
point correlation function, or its Fourier transform, the bis- 
pectrum are often used. The higher-order statistics, however, 
are only part of the whole picture about non-Gaussianity. It 
takes a full hierarchy of n-point correlation functions, or the 
polyspectra, to complete the statistical characterisation of 
the CMB anisotropies. 

Since the release of 1-year WMAP data, great ef- 
forts have been made for the search and detection of non- 
Gaussianity via various approaches and methods (Chiang et 
al. 2003; Park 2004; Eriksen et al. 2004c; Vielva et al. 2004; 
Copi, Huterer & Starkman 2004; Eriksen et al. 2004b; Ca- 
bella et al. 2004; Hansen et al. 2004; Mukherjee & Wang 
2004; Larson & Wandelt 2004). These detection of non- 
Gaussianity shall augment the error bars on the estimation 
in cosmological parameters. 

Practically speaking, characterizing non-Gaussianity 
through phases is one of the most general approaches. Based 
on the random-phase hypothesis, the key to developing sta- 
tistical methods using phases is testing their "randomness" . 
Due to the 2n wrapping, however, the phases tend to be uni- 
formly distributed at [0, 2n]. For example, for a point source 
produces phases are distributed evenly and orderly, but not 
randomly between and 2-k. Probing non-Gaussianity via 
examining uniformity of phases themselves between and 2ir 
is often as ineffective as via examining one-point Gaussian 
(temperature) probability distribution p(T) (e.g. it is possi- 
ble p(T) is still Gaussian for a non-Gaussian field). The 2ir 
wrapping causing a uniform distribution of phases is similar 
to the central limit theorem in action producing one-point 
Gaussian probability distribution. We thus seek associations 
between phases as a more sensitive and effective statistical 
measure. 

The linear association such as auto correlation (</></> ), 
however, does not give useful statistics due to the circular 
nature of phases. To counter this problem, return mapping of 
phases is introduced to associate phase pairs systematically 
(Chiang, Coles & Naselsky 2002). The main idea is mapping 
all phase pairs with the same separation (A£, Am) onto a 
square, which is conceptually similar to the auto correlation 
£(A£, Am) = {4> e , m <t>e+&i, m +&m)- 

Another important feature of phases is that phases are 
closely related to morphology. Through pixel-by-pixel cross 



correlation, maps with the same phases display strong re- 
semblance in morphology, regardless of their power spec- 
trum (Chiang 2001). The level of cross-correlation of phases 
taken from two images therefore renders significance of re- 
semblance between them. Based on the simple but prevail- 
ing assumption that the CMB signals should not correlate 
with the foregrounds (i.e. the microwave foregrounds should 
not have knowledge in what the CMB signals 'look like'), 
cross correlations of phases between the derived CMB and 
the foregrounds shed light on the status of microwave fore- 
ground cleaning. Dineen and Coles (2003) perform cross cor- 
relations in pixel domain between the derived CMB and the 
foreground maps. Naselsky, Doroshkevich & Verkhodanov 
(2003; 2004), Naselsky et al.(2004) use cross correlation of 
phases to illustrate the foreground contaminations in the 
derived CMB signals. 

Return mapping of phases renders phase associations on 
a square with each side ranged [0, 2n] (with periodic bound- 
aries) . Cross correlation of phases between derived CMB and 
foregrounds connect phases from the same aim mode of the 
maps also produce mapping in a [0, 2-k] square. The null 
hypotheses for both cases: random phases for a Gaussian 
CMB sky and no cross correlation between the CMB and the 
foregrounds shall result in random points on [0, 2ir] squares. 
Statistics that are developed for a square taking into account 
of the periodic boundaries can be implemented both on the 
return mapping (auto correlation of phases) of the derived 
CMB signal and on cross correlation of phases between CMB 
and foregrounds. More importantly, through checking both 
auto and cross correlations of phases from the derived maps, 
one can gain insight into the issue of non-Gaussianity. 

In this paper we test on each spherical harmonic multi- 
pole t the auto and cross correlations, as well as the overall 
correlations. This paper is arranged as follows. In Section 
2 we recap the phase mapping technique and introduce the 
Kuiper's statistics. We apply the Kuiper's statistics in Sec- 
tion 3 both on the auto and cross correlation of phases of 
the 1-year WMAP derived maps and connect the phase cor- 
relations with non-trivial whitened AT distribution in ID 
Fourier composition. The conclusion and discussions are in 
Section 4. 



2 THE PHASE MAPPING TECHNIQUE AND 
THE KUIPER'S STATISTICS 

2.1 The return mapping of phases and 
two-dimensional uniformity 

Chiang, Coles & Naselsky (2002) have introduced a tech- 
nique called return mapping of phases to render associa- 
tions between phase pairs on a square. The idea is borrowed 
from the return map in chaotic dynamics (May 1976). The 
phase pairs are taken systematically and mapped onto a 
[0, 2n] square. In a return map of phase pairs with the sep- 
aration (A£, Am), for example, the mapped points formed 
from phase pairs have the coordinates (cf>e,m , 0<+Af,m+Am) 
for all possible i and m. Because of the circular nature of 
phases, the [0, 2n] return map is periodic at all 4 sides, which 
can be viewed as a flat torus. 

After the mapping of phases, various statistics can be 
applied to extract the statistical significance. Chiang, Nasel- 
sky & Coles (2004) use a simple mean chi-square statistic to 
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Figure 1. Wc show two panels formed with a mosaic of 4 return 
maps of phases (left) and the projection of the mapping points to 
tackle the circular nature of phases. On the left top, the 4 points 
(3, 7, e and i) on the return maps in fact represent the same point 
due to periodicity of the return maps. To probe the connectivity 
of points on the return maps, therefore, we consider the simplest 
directions: anti-diagonal (BD), and diagonal (AC) projecting onto 
the y axis. On the right panel, we put 3 points a, b and c, as an 
example on the return map and probe the connectivity in the 
anti-diagonal direction. We reproduce the same map on top of it: 
a , 6 and c . Due to the projection, we consider the same total 
area but different configuration (the shaded area), which covers 
the 2 points from the lower map (a and 6), and 1 from the upper 
map (c ). The 3 seemingly non-correlated points in a single return 
map now show the alignment along the anti-diagonal direction. 



extract the information of the one- dimensional uniformity 
of the distribution of the return maps. In order to probe 
associations between points (uniformity in 2 dimensions), 
care has to be taken on the periodicity of the return maps. 
In Fig.l on the left we show the mosaic of 4 return maps to 
indicate the periodicity. To probe the connectivity of points 
on the return maps, we consider the following 2 directions: 
anti-diagonal (BD), and diagonal (AC) onto the y axis. This 
projection is also useful for cross correlation of phases. Com- 
plete cross correlation of phases produces points exactly on 
the anti-diagonal line. 

The projection of a point located at (xo,yo) along the 
anti-diagonal direction to y axis as shown in Fig.l is equiv- 
alent to taking the difference of its x and y coordinates, 
i.e. yo — xo; in the diagonal direction yo + xo- In a re- 
turn map of phases with separation (Af, Am), points are 
formed with phase pairs (f>i m and 4't+Ai.m+Am, the projec- 
tion then becomes (j>e+Ae, m +Am — <t>e m '- the phase difference. 
In Coles et al. (2003) they probe phase correlations by inves- 
tigating the uniformity of the neighbouring phase difference 
4>i,m+i — 4>im- This is equivalent to testing the uniformity 
of the projected points on a return map of the separation 
(A£, Am) = (0, 1). Following this line of thought, the anal- 
ysis of the projection in either the vertical or the horizontal 
directions is that of the randomness of <j>e m themselves. 
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Figure 2. The difference in the statistical significance of Kuipcr 
statistics from its analytical derivation, Eq.(5) and (6), and Monte 
Carlo simulation. Wc list the 3 frequently used significance value 
(from top to bottom): 66.269%, 95.450% and 99.730%, respec- 
tively, which corresponds to Gaussian 1, 2, 3-cr C.L.. The y-axis 
is the difference in absolute value, and the x-axis is the data point 
N. 



2.2 The Kuiper's statistic 

The Kuiper's statistic can be viewed as a variant of 
Kolmogorov-Smirnov (K-S) test (Kuiper 1960; Press et al. 
1992; Fisher 1993). As the projection of points on the di- 
rections we mention in §2.1 produces unbinned distribution 
that is a function of single independent variable, it is very 
useful to apply the K-S test to probe its uniformity. The K-S 
statistic is taken as the maximum distance of the cumulative 
proability distribution against the theoretical one: 



D 



max \Sn(x) 

-oo<a:<oo 



P(X)\ 



(3) 



For circular function, however, one needs to take into ac- 
count of the maximum distance both above and below the 
P(x) 



V 



D+ + D- = 
+ 



max [Sn(x) 

-oo<txKoo 

max [P{x) - 

-O0<3?<00 



" P(x)] 
S N (x)] 



(4) 



and the C.L. against the null hypothesis (e.g. uniformity of 
phases in our case) can be calculated from (Press et al. 1992) 



VN + 0.155 + ^= 



C.L. = 1 — QKuipcr yV 
where 

oo 

Q Kui pc r (A) = 2^(4/ A 2 - l)e- 2j2A2 
i=i 

and N is the data points. 

We run Monte Carlo (MC) simulation in order to test 
Eq.(5) and (6). One hundred thousand realizations of the dif- 



(5) 



(6) 
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map (WFM) 3 by Tegmark, de Oliveira-Costa & Hamilton 
(2003), and the Phase-Cleaned Map (PCM) by Naselsky et 
al. 4 (2003). Note that these maps are used only for testing 
the effectiveness of the Kuiper's statistics because of their 
different approaches on foreground cleaning and their dif- 
ferent morphology. Our claim about non-Gaussianity in this 
paper shall not contradict with others in the literature as 
most of these maps, as advised by the authors, are not for 
scientific purposes. 

In Fig. 3 as an example of how the Kuiper's statistics 
works, we show the cross correlation of phases between the 
ILC and the W-band foreground map at £ = 3 (top) and 
£ — 8 (bottom) and their corresponding cumulative proba- 
bility functions from the projection in the anti-diagonal di- 
rection onto the y axis. We show with thin (D+) and thick 
(D-) shaded lines the maximum distances of the cumulative 
probability distribution above and below the theoretical one, 
respectively. 



12 3 4 5 
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Figure 3. The cross correlation of phases at t = 3 (top-left: 
m = — 3) and 8 (bottom-left: m = — 8) between the ILC 
map and the W-band foreground map. One can imagine that, for 
example, for t = 3 the 4 points are highly correlated along the 
anti-diagonal direction if the left-most point is repositioned at 
the far-right. The same for £ = 8 if we reposition the 3 points at 
the lower-right region to the above top axis line (cf Fig.l). The 
panels on the right are their corresponding cumulative probability 
functions. The dashed line is the P(x), the theoretical cumulative 
probability distribution for 2D uniformity. The thin and thick 
shaded lines are the D+ and D_ , respectively, their maximum 
distances of the cumulative probability distribution above and 
below the theoretical one. 



ference of two random phase series (the phases are defined 
[0, 2n)) are simulated, which is to mimic the projection of 
return map shown in Fig.l. For each N, we have an ensemble 
of 100 000 random realization each producing maximum dis- 
tance V as in Eq.(4). We then sort these 100 000 V to yield 
the statistical significance from MC simulation. Simultane- 
ously, the statistical significance can be calculated directly 
from Eq.(5) and (6). In Fig. 2 we show the difference in the 
statistical significance of the Kuiper statistics from its ana- 
lytical derivation, Eq.(5) and (6), and MC simulation. We 
list the 3 frequently quoted significance value (from top to 
bottom): 66.269%, 95.450% and 99.730%, which corresponds 
to Gaussian 1, 2, 3-ct C.L., respectively. One can see that for 
95.450% and 99.730% the difference between analytical cal- 
culation and MC simulation is down to 10~ 3 level even for 
as few data points as N = 10. 

We use the maps available on the WWW for testing 
this method: the derived foreground maps at 5 frequency 
bands from the WMAP website 1 and the 4 derived CMB 
maps : the internal linear combination (ILC) map by the 
WMAP science team (Bennett et al. 2003a; Bennett et al. 
2003b; Bennett et al. 2003c; Hinshaw et al. 2003a; Hin- 
shaw et al. 2003b; Komatsu et al. 2003), the ILC map by 
Eriksen et al. (hereafter EILC) (2004a) 2 , the Wiener-filtered 



3 AUTO AND CROSS CORRELATION OF 
PHASES 

3.1 Auto correlation 

We consider first within each multipole number £ the map- 
ping of phases of fixed Am. In Fig. 4 we show the Kuiper's 
statistics on return maps for the WMAP ILC map, the EILC 
map, the WFM and the PCM. Within each multipole £ we 
take the phases of fixed separation Am and map them into 
a return map. For the multipole number I — £o we take the 
separation from Am = 1 to Am = £o/2, although in general 
we can take any Am values to test the randomness of phases. 
In order to display the C.L. more clearly against random 
phase hypothesis, we show the QKuipcr instead. So 10~ 3 cor- 
responds to C.L. 99.9% against uniformity of phases, 10 -2 
corresponds to C.L. 99% . . . etc.. 

One interesting result is that for £ — 6 the mapping of 
Am = 1 — 3 for all 4 maps are all above 75% C.L.. and the 
multipole numbers which has Am mappings that are above 
95.45% C.L. against random phase hypothesis are I = 6, 12, 
15, 17, 19, 30, 35 - 36, 38 - 39, 41, 43, 45 - 47, 49 - 50. For 
Am = 6 at t = 38 and Am = 5 at £ = 46, both mappings 
reach above 99.73%. For other 3 maps, they all have modes 
that reach 95.45%. We do not, however, claim the overall 
non-Gaussianity of these maps. 

We also consider the mapping for separation (A£, Am) 
from all the available phases from £ = 2 to 50 (except those 
of m = modes). In Fig. 5 we show the Kuiper's statistics 
for the 4 maps: the ILC, EILC, WFM and PCM. The sep- 
aration of return mapping ranges from to 5 for both A£ 
(x axis), Am y axis), although it can be extended to higher 
number. The ILC, EILC and PCM have (A£, Am) = (2,0) 
correlation above 68.27% C.L. and there are considerable 
degrees of phase correlation at separation ( Al, 0) for ILC, 
EILC and PCM, but except for WFM. 

Note also the strong correlation of (A£, Am) = (4, 0) for 
the PCM. This specific correlation is recently mentioned in 



1 http://lambda.gsfc.nasa.gov 

2 http : //www . astro . uio . no/~hke/cmbdata/ 



3 http : //www . hep . upenn . edu/~max/wmap . html 

4 http : //www . nbl . dk/~chlang/wmap . html 
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Figure 4. The Kuiper's statistics on return maps for the WMAP 
ILC map, the EILC map, TOH's Wiener-filter map and the PCM 
map. In order to display the C.L. more clearly against random 
phase hypothesis, we show the QKuiper instead. So 10~ 3 corre- 
sponds to C.L. 99.9% against uniformity of phases, 10~ 2 corre- 
sponds to C.L. 99% . . . etc.. Within each multipole I we take the 
phases of fixed separation Am and map them in a return map. 
We take Am = 1 up to Am = lg/2 for multipole number Iq. 
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Figure 5. The Kuiper's statistics on return map for separation 
(Af, Am). The mapping takes all the phases from I = 2 to 50. 
The separation of mapping is ranged to 5 for both 81 and 8m. 
In order to display more clearly the C.L. against random phase 
hypothesis, we show the QKuiper instead. So lO -3 corresponds 
to 99.9% C.L. against uniformity of phases, 10~ 2 to 99% C.L. 
. . . etc.. 



Naselsky & Novikov (2005), which is caused by the symmet- 
ric (w.r.t. Galactic centre) point-like peaks lying along the 
Galactic plane. The PCM is produced without further clean- 
ing performed for ILC, EILC and WFM, hence the Galactic 
contamination is still present. 

In order to see how non-random phases affect the dis- 
tribution, we assemble only the phase part in each harmonic 
number £ in Fourier composition (the so-called "whitened" 



AT* , 



E 



at n 



■ exp(imx) 



exp(ii?W) exp(ima;). 



(7) 



In Fig. 6 we display for I = 2 to 17 the "whitened" distri- 
butions which have no influence from the amplitudes |o^ m |. 
This representation in ID Fourier composition for each £, 
rather than in the standard whole-sky spherical harmonic 
composition, can provide some advantages for further anal- 
ysis. The morphology can be seen more clearly in ID Fourier 
composition than in spherical harmonic composition. Here 
we can relate phase correlation (in Fourier space) to the sim- 
ple GRF peak statistics (in real space). The extrema (max- 
imum or minimum) that have reached 3u appear in the dis- 
tribution of I — 6, 12, 15 and 17: —3.0a, 3.1ct, — 3.2cr and 
3.3cr, respectively. These multipole numbers correspond to 
those picked up by the Kuiper's statistics with non-random 
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Figure 6. The assembled "whitened" distribution I by i from the WMAP ILC map. We assemble only the phase part exp(i<f)£ m ) of each 
harmonic number £ in Fourier composition: AT( , whitened C^) = y"] m exp(i<f>£ m ) exp(imx). The extrema that reach 3<r are i = 6, 12, 15 
and 17: —3. Oct, 3.1ct, — 3.2a and 3.3<r, respectively. These multiple numbers are picked up by phase correlations with C.L. 95.45% from 
the Kuiper's statistics. 



phases at 95.45% C.L.. Furthermore, phase correlation does 
not only manifest itself in the extrema. One can see, for 
example, for I = 10 the assembled whitened AT has re- 
peated peaks. The C.L. against random phase hypothesis at 
Am = 2 of £ = 10 reaches 94.13%. 

3.2 Cross correlation 

The Kuiper's statistics can also be applied to cross correla- 
tion. 

In Fig. 7 we show the C.L. against the null hypothesis in 
cross-correlation of phases for £ = 2 to 50. The thick lines are 
the WMAP ILC cross the WMAP foreground maps at (from 
top to bottom) K, Ka, Q, V and W-band, the cross sign (x), 
the diamond sign (o), the plus sign (+) are the EILC by 
Eriksen et al. (2004a), the PCM by Naselsky et al.(2003) and 
the WFM by Tegmark et al.(2003), respectively, cross the 
WMAP foreground maps. We would like to point out that 
the ILC are produced by less weights from the W-band map 



than V and Q. One can see that the ILC map has substantial 
cross correlation of phases for I = 3 and 8 with WMAP all 
5 channel foreground maps. High cross correlation with the 
foregrounds, however, does not necessarily imply the non- 
Gaussianity for that multipole number. 

In Fig. 8 we show the assembled AT distribution for £ — 
8 from WMAP ILC map and WMAP W-channel foreground 
map. From the Kuiper's statistics for cross correlation of 
phases, for I = 8 the correlation reaches C.L. 89.98%, which 
reflects on their resemblance in morphology. 



4 CONCLUSION AND DISCUSSIONS 

In this paper we introduce the Kuiper's statistic to probe 
the 2D uniformity for the phase mapping technique. The 
Kuiper's statistics are useful to explore both auto correla- 
tion and cross correlation of phases. We use the 4 maps to 
test the effectiveness of this method: WMAP ILC, the ILC 
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by Eriksen et al., the WFM by Tegmark et aland the PCM 
by Naselsky et al.and found several multipole numbers with 
non-randomness of phases over 95.45%. The Kuiper's statis- 
tics can also used to test on all the available phases. 

Contrary to the AT representation for each I from stan- 
dard spherical harmonic composition, we use ID Fourier 
composition to display the AT. We use the "whitened" 
Fourier composition to display the connection between non- 
random phases and non-trivial morphology, and between 
cross-correlated phases of 2 maps and resemblance in their 
morphology. 

The peculiarity of £ = 3 and 8 found by Copi et al.(2004) 
can be seen clearly in the cross correlation with the fore- 
grounds at all 5 foreground maps, as shown in Fig. 7, par- 
ticularly with the W band foreground map. This is an- 
other advantage of using phases as cross checking for non- 
Gaussianity for the CMB signal. There is no other method 
so far that can cross check both the foreground maps and 
the derived CMB map when peculiarities are found. What 
is unclear though, is the non-random phases at I = 6 ap- 
pearing on all ILC, EILC WFM, and PCM (Fig.4), where 
the cross correlation with foregrounds is reasonably low. 
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